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We study, in the dynamically evolving QGP fireball formed in relativistic heavy ion collisions 
at RHIC and LHC, the growth of strangeness yield toward and beyond the chemical equilibrium. 
We account for the contribution of the direct strangeness production and evaluate the thermal- 
QCD strangeness production mechanisms. The specific yield of strangeness per entropy, s/5 1 , is the 
primary target variable. We explore the effect of collision impact parameter, i.e., fireball size, on 
kinetic strangeness chemical equilibration in QGP. Insights gained in study the RHIC data with 
regard to the dynamics of the fireball are applied to the study strangeness production at the LHC. 
We use these results and consider the strange hadron relative particle yields at RHIC and LHC in 
a systematic fashion. We consider both the dependence on s/S and directly participant number 
dependence. 



PACS numbers: 25.75.Nq, 25.75.-q, 24.10.Pa,12.38.Mh 

I. INTRODUCTION 

Conversion of kinetic collision energy into high mul- 
tiplicity of newly made hadronic particles is one of the 
most notable features of reactions observed at the Rel- 
ativistic Heavy Ion Collider (RHIC) at the Brookhavcn 
National Laboratory (BNL) [1]. In this process, aside 
of the light u and d quark pairs, present in all matter 
surrounding us, the strange flavor quark pairs s, s are 
produced copiously (in general, in what follows, the par- 
ticle symbol will refer to the corresponding particle yield, 
either total or per unit of rapidity, as appropriate). The 
final s-yield depends on the initial reactions, and on the 
history of the fireball, and thus, also on the nature and 
properties of the phase of matter formed. On the time 
scale of hadronic interactions, strangeness flavor is con- 
served, and prior to any weak interaction decays, we have 
s = s — when we refer to strangeness yield, production, 
etc, we always address yield, production, etc, of strange 
quark pairs. 

This study is addressing strangeness under the physical 
conditions achieved, at RHIC, at the highest attainable 
reaction energy today, ^/snn = 200 GeV, and in future 
at LHC. We are particularly interested in the sensitiv- 
ity of strangeness production to the nature and prop- 
erties of the matter formed in the heavy ion reactions. 
Theoretical studies have shown that strangeness is pro- 
duced rapidly in collisions (fusion) of thermalized glu- 
ons @, [H, within the deconfmed state, the quark-gluon 
plasma (QGP) formed in the central collisions of heavy 
nuclei. On the scale of RHIC reaction time r < 10 fm, 
the hadron based reactions were found to be too slow to 
allow copious strangeness production after thermaliza- 
tion of matter and are even more ineffective to produce 
multi-strange hadrons ||. 

On the other hand, strangeness can be produced fast 



in the QGP phase, as we shall see achieving near chem- 
ical equilibrium in QGP phase formed at RHIC and 
even overshooting the chemical equilibrium at hadroniza- 
tion at the LHC. Thus, the situation is quite differ- 
ent when the deconfmed QGP state breaks up in a 
fast hadronization process: the enhancement of strange 
hadrons and most specifically strange antibaryons, grow- 
ing with valance strange quark content of hadrons pro- 
duced is the predicted characteristic property of the de- 
confined QGP phase This happens since in the 
breakup of the strangeness rich deconfmed state, i.e., 
hadronization, several strange quarks formed in prior, 
and independent, reactions can combine into a multi- 
strange hadron. 

Our main objective in this work is to quantify the 
mechanisms of kinetic strangeness production occurring 
in thermal gluon collision (fusion) processes during the 
expansion phase of the quark-gluon fireball. We ex- 
plore the centrality dependence in a wide range between 
peripheral, and most central reactions, in which up to 
90% of projectile and target nucleons participate. In 
this regard, this paper is a theoretical companion to the 
more phenomenological analysis of experimental RHIC 
data 0, HI, and uses the insights gained in this analysis, 
in particular, regarding the dynamics of the QGP expan- 
sion. This is then applied to extrapolate our approach to 
the LHC energy domain, where our prior particle yield 
study was based on a parametric consideration of final 
state strangeness yield Q. 

We also study in depth the centrality dependence of 
strangeness production in the RHIC-LHC energy range. 
As the centrality of the nuclear reaction and the number 
of participants A decreases, the number of thermal col- 
lisions gradually diminishes, and with it the strangeness 
enhancement effect also diminishes gradually, similarly to 
the AGS-SPS energy range [1,01 • This decrease drives in 
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turn a gradual decrease in the centrality dependent pro- 
duction rate of multistrange hadrons. This behavior we 
discuss in detail here is an important and characteristic 
phenomenological feature of the kinetic particle collision 
mechanism of strangeness production and enhancement. 

Interestingly, in this regard the kinetic mechanism 
of strangeness production differs from models, such as 
'canonical enhancement model', which are deriving the 
strange hadron enhancement as a result of an always 
prevailing hadronic phase chemical equilibrium 

[13, Hlll- 
Tlie volume dependence of the canonical phase space 
yields [l2j, and the smallness of the N-N reference sys - 
tems produce the centrality enhancement effect |13j . 
However, 'canonical enhancement' is rising very rapidly 
considering rather small collision systems, and with de- 
creasing energy [HI ]. 

The understanding of strangeness production during 
the expansion phase of the quark-gluon fireball allows 
us to study in depth the reaction mechanisms which are 
determining the final state yield of strangeness. In this 
way, we learn how 'deep' into the history of QGP expan- 
sion this observable allows us to look. Collective matter 
flow features observed at RIIIC suggest that thermaliza- 
tion of parton matter occurred very fast, i.e., the en- 
tropy S has been produced in a not yet fully understood 
fast process of parton thermalization, prior to the pro- 
duction of strangeness pair yield s. We formulate the 
kinetic equations allowing to address, in some detail, the 
growth in specific strangeness per (fixed) entropy s/S in 
the thermal QGP processes in both a longitudinally and 
transversely expanding QGP fireball. 

The time evolution of s/S has been considered for 
the first time early on in the development of the QGP 
physics [HI]. However, the model of dense matter evolu- 
tion, and the range of statistical parameters considered 
at the time is not appropriate for the RHIC and LHC 
physics environments. Moreover, we recompute here the 
rate of strangeness production knowing the best current 
values of QCD parameters, the coupling constant a s , and 
the strange quark mass m s . These parameters alter de- 
cisively the values of s/S in chemical equilibrium, and 
thus the dynamical time scale of the approach to chemi- 
cal equilibrium. In order to be able to compute the evo- 
lution in time of s/S, we must also evaluate how near to 
chemical equilibrium is strangeness in QGP: this near- 
ness is characterized by a parameter 7 S , roughly the ratio 
of prevailing strangeness density to chemical equilibrium 
density at prevailing temperature. 

Setting up the production of strangeness, we assume 
here that it follows in time the chemical equilibration of 
the light q — u, d quarks, and g gluons [la |. which are 
believed to occur at 1 fm scale. These effectively massless 
particles can be produced by entirely soft processes which 
are intrinsically non-perturbative [171 ]. Their chemical 
equilibration can be further driven by mult i- particle col- 
lisions [H, [l!|. Considering these studies, we assume 
here relatively short relaxation times for q, g, we can- 
not compute these using the same perturbative method 



as will be developed here to evaluate strangeness yield 
equilibration. It is the finite strangeness mass combined 
with the measured strength of the running QCD coupling 
constant a s which allows us the use of perturbative for- 
malism in study of strangeness production [20| with some 
minimal confidence. 

The chemical relaxation times for the strangeness ap- 
proach to chemical equilibrium, in an expandin g Q GP, 
has been considered several times before [!, Hi], l22l l23l 
HH HH, Hil • The study of QGP strangeness chemical equi- 
libration must not be confused with the phenomenologi- 
cal investigation of chemical equilibrium in the final state 
hadron abundance. Because hadron phase space is gener- 
ally smaller, chemical equilibrium and indeed excess over 
equilibrium is much more easy to attain, and there is 
some continuing discussion of this question [U [H, [29[ . 

We offer in our work a comprehensive exploration how 
the impact parameter dependence, and consideration of 
energy dependence, influences chemical equilibration in 
QGP. We evolve in time not the strangeness itself, but 
the specific strangeness per entropy s/S. In this way, 
we can identify more clearly the production processes of 
strangeness, since entropy is produced earlier on, and 
is (nearly) conserved during the time period of thermal 
strangeness production. 

Importantly, s/S is an experimental observable, prac- 
tically preserved in the fast hadronization process. Thus, 
we can connect the final state of the QGP evolution di- 
rectly to experimental soft hadron yield experimental re- 
sults. We find that much of the variability about the 
initial conditions, such as dependence on initial temper- 
ature, cancels in the final result. This specific observable 
will be shown to yield nearly model independent insights 
about thermal strangeness production in QGP. 

In the following section |TT1 we begin with a brief discus- 
sion of general features relevant in all considerations pre- 
sented. In subsection III D[ we formulate the kinetic equa- 
tions describing the growth in specific strangeness per 
entropy s/S and show that, at RHIC, the observed spe- 
cific per entropy strangeness yield suggests that the 'di- 
rect' and 'thermal' processes are of comparable strength 
in most central reactions. We discuss the magnitude 
and importance of QCD parameters in subsection III El 
In order to integrate as function of reaction time the 
strangeness yield equations, we develop a simple collec- 
tive expansion model of the plasma phase in subsection 

Em 

In section IIIII we study the thermal strangeness pro- 
duction processes. We then obtain reference yields for 
two different expansion geometries at RHIC in subsec- 
tion [TTI3 We extrapolate this to the LHC environment 
in subsection IIII Bl We explore how 'deep' into the early 
history of the hot and dense fireball the strangeness sig- 
nature of thermal QGP is allowing us to look, i.e., the 
dependence on initial conditions, in subsection llll CI - im- 
portantly we find that the selection of the initial value of 
s/S related to direct production of strangeness, has only 
minor impact on the final results regarding strangeness 
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yield. We then explore the influence of fundamental un- 
certainties, such as the present day limited knowledge 
about the strange quark mass and the QCD thermal ef- 
fects on the freezing of the strange quark degrees of free- 
dom, in subsection IIII Dl 

In section IIV1 we connect the results we obtained to 
the experimental particle yields. In subsection IIV Al 
we consider the relationship between hadron multiplic- 
ity and entropy yield, and obtain the strange hadron 
yield as function of s/S. This allows an assessment 
how strangeness yield, at RHIC and LHC, influence 
physical observables. In particular, we discuss how 
K+/7T+ changes between these two experimental en- 
vironments. Then, we discuss, in subsection IIVBI 
for the two most often used statistical hadronization 
models (sudden hadronization and chemical equilibrium 
hadronic gas (HG) hadronization) the production of 
strange hadrons as function of participant number A, 
keeping the hadronization condition independent of A. 
In subsection IIV C[ we apply the insights gained in study 
of thermal strangeness production to evaluate thermal 
charm production at the RHIC and the LHC environ- 
ments. 



II. REMARKS ABOUT STRANGENESS 
PRODUCTION AND DENSITY 

A. Parton equilibration and strangeness 
production 

The total final state hadron multiplicity is a measure 
of the entropy S produced prior to thermal production 
of strangeness s in QGP: once a quasi-thermal exponen- 
tial energy distribution of partons has been formed, the 
entropy production has been completed. Further evo- 
lution of the dense deconfined fireball is nearly entropy 
conserving, even though it is strangeness flavor produc- 
ing: fusion of gluons, or light quark pair annihilation into 
strangeness, is a nearly entropy conserving process [30j . 
We note that, in reactions between two thermal particles 
into a strangeness pair, the energy content of each initial 
state parton is transferred to the two reaction products, 
so thermal partons produce thermal shape of strangeness 
spectrum. 

The entropy produced, in RHIC reactions, has been 
evaluated in recent studies of hadron multiplicities. In 
Au-Au reactions, at ^/snn = 200 GeV, one sees S ~ 
35,000. Furthermore, at central rapidity, the yield of 
entropy is dS/dy ~ 5000. In the benchmark results we 
present for LHC, we will assume that the the central 
rapidity entropy yield is about 4 times greater than at 
RHIC. 

The temporal evolution of the QGP fireball ends when 
the temperature has decreased to the QGP hadroniza- 
tion value. In the breakup of the QGP, the yields of 
hadrons are established, and it is rather difficult in the 
ensuing rather short lived evolution lasting not more 



than 1.5 fm/c to alter these yields appreciably. Thus, 
the hadronization volume, with the normalizing factor 
dV/dy, provides the normalization of hadron particle 
yields per unit of rapidity. The final value of dVirt)/ dy is 
result of analysis of hadron particle yields and our model 
of the time dependence of dV{r)/dy will be constrained 
by the magnitude of dVf / dy obtained in Ref . Q • 

There are two separate stages of strangeness (charm) 
production, corresponding to the two practically distinct 
periods of the fireball evolution: 

a) 'direct' production creates a 'background' yield corre- 
sponding to what might be obtained in a superposition 
model of independent nucleon-nucleon (N-N) reactions. 

b) strangeness production in collisions between thermally 
equilibrated QGP constituents. 

For b) to be relevant prior formation of a thermal, de- 
confined QGP phase is required. Without process b), 
the yield of strangeness should not be enhanced in A-A 
collisions as compared to scaled N-N reactions. We will 
not discuss in detail mechanisms a) of direct particle pro- 
duction here, our interest is restricted to (approximate) 
initial yields, which are the baseline for the thermal mech- 
anisms acting in the QGP, and define the strength of any 
enhancement. 

We note for the record, that strangeness initial pro- 
duction, like other relatively soft parton production pro- 
cesses, are believed to be due to color string breaking 
mechanism [3l| and the Pythia 6 model of soft hadron 
production presumes that the relative strength of u : d : s 
production is 1 : 1 : 0.3. The (initial) charm production 
is due to high energy parton collisions [321 ]. 

Our study of the thermal strange particle production 
processes is based on kinetic theory of particle collisions. 
There is considerable uncertainty about the initial mo- 
mentum distributions of soft partons present in the initial 
state. However, two recent theoretical studies argue that 
there is rapid thermalization. The nonlinear gluon pro- 
duction processes leads to the gluon momentum distri- 
bution equilibration [19]. Axial asymmetry of the initial 
state causes collective instabilities which further acceler- 
ate thermalization of partons [331 ]. 

In this context, it is important to advance one result of 
our study, namely that there is little sensitivity to the ini- 
tial thermal condition: a wide range of 'reasonable' initial 
temperatures leading to very similar strangeness produc- 
tion results, as long as the entropy content is preserved. 
In order to understand this, consider a decrease in initial 
temperature. This requires, at fixed entropy, an increase 
in initial volume, and this, then, is associated with in- 
creased lifespan of the fireball in the high temperature 
strangeness producing domain. These two effects com- 
bine to compensate the reduced strangeness production 
rate per unit of time and volume that is associated with 
reduced ambient temperature. 

We believe that this mechanism also implies that the 
precise form of the momentum spectrum of the initial 
state partons is of minor practical relevance for the pur- 
pose of evaluation of strangeness production, and we do 
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not study this. Therefore, without loss of generality, we 
can assume that the parton distributions we use in the ki- 
netic strangeness formation process have thermal shape, 
and the ambient temperature is determined considering 
(lattice fitted) equations of state relation of initial tem- 
perature and entropy density [13], for a given geometric 
initial volume. 

We already remarked above that the production of 
strangeness, in a cascade of N-N reactions (without de- 
confinement), is not able to add significantly to the initial 
strangeness yield considering the short lifespan of the fire- 
ball. Thus, this alternative will not be further considered 
in this work. Similarly, any additional strangeness pro- 
duced in the rapid hadronization of QGP into hadrons 
must be negligible compared to the thermal production 
process which occurs at higher particle density (temper- 
ature), and during a considerably longer lifespan. 

B. Approach to chemical yield equilibrium 

It has been shown, considering the entropy maximiza- 
tion principle, that the approach of particle densities to 
chemical equilibrium density can be characterized by 
the statistical parameter 7* [35j, which varies with the 
local proper time r during the collision. For example, for 
gluons, 

f ( T \ p -E/T 

Pg (r)^J d 3 P l l^ (r)e _ E/T , E=^T?. (1) 

Generally, the Lagrangian mass of gluons is zero. How- 
ever, one may be tempted to think that thermal mass 
m(T) could change decisively results, suppressing the col- 
lisional strangeness production. However, one finds that 
instead the process of gluon decay becomes relevant and 
if at all, there is a net rate increase of strangeness pro- 
duction [36[. One could argue that the scheme to study 
kinetic process of strangeness chemical equilibration us- 
ing thermal mass amounts to a different resummation of 
reaction processes. In this work, we will consider the evo- 
lution of 7^ gp (t) based on Lagrangian masses, allowing 
for 7°- gp (t) and 7^ GP (r). 

For strange quarks, we will keep only the Boltzmann 
term (ignoring the denominator in Eq. ([I]) above) and 
thus: 

Ps (r) 7 ? GP (r) P r = 7? GP f^(*), z = ^. (2) 

Here, g s is the strange quark degeneracy, and K2 is a 
Bessel function. Since we will employ strangeness occu- 
pancy in hadron phase, 7^, we have included the sub- 
script QGP. We will henceforth drop this subscript, and 
occupancy parameters without superscript will, in gen- 
eral, refer to the QGP phase, while the hadronic gas 
phase variable, when these are expected to differ from 
QGP, will have a superscript h. 

Our target variable is the final QGP state specific yield 
of strangeness per entropy, s/S, and the related phase 



space occupancy j s . Both these variables have an impor- 
tant physical relevance: s/S determines the final yield of 
strange hadrons compared to all hadrons, and its value 
implies some particular yield of reference yields, such as, 
e.g., K + /7r + . 7 S characterizes the approach to chemi- 
cal equilibrium, it measures strangeness yield in terms of 
the chemical equilibrium yield. The strangeness phase 
space of QGP and HG phases are different. Strangeness 
in QGP is much denser than in the HG phase, con- 
sidering the range of strange quark masses, 0.080 < 
m s (fi = 2GeV) < 0.125 GeV. Therefore, (near) chem- 
ical strangeness equilibrium in the QGP phase 7 S ~ 1, 
implies a significantly oversaturated hadron phase space 
abundance y~ > 1 after hadronization. 7^ is directly 
controlling the relative yields of hadrons with different 
s + s valance quark content and is thus observable. 



C. Role of initial conditions 

We do not understand well the conditions in the QGP 
phase at time as early as To = 0.25 fm for RHIC, and 
To = 0.1 fm for LHC, when we presume that the thermal 
momentum distribution is practically established. Thus, 
we must make a number of assumptions and check if these 
impact our results. The relevant parameters that could 
govern the strangeness production are: 
7g(7") and, in particular, the initial value at To; 
To, the time at which we assume thermal momentum of 
partons is reached; 

s/S\ To is the initial strangeness yield originating in direct 
parton collisions; 

R± is the transverse radius dimension at initial time, re- 
lated to the collision geometry; 

(t) is the transverse expansion velocity, and in partic- 
ular its maximum value at hadronization; 
T g and T q , the relaxation time constant of gluon and 
quark fugacities, considering that quarks are less relevant 
compared to gluons with regard to strangeness produc- 
tion, we will assume t 9 = 1.5 t s throughout this work. 

We explore, here, a characteristic gluon thermaliza- 
tion time 0.1 < To < 1.5 fm/c, with the longest value 
applicable to most peripheral RHIC Au-Au reactions at 
^/snn = 200 GeV, and the shortest period assumed for 
the future LHC central collisions. Knowing the exact dy- 
namics of thermalization and how long it takes will be, 
as we shall see, rather unimportant for the final insights 
we obtain. 

All initial state parameters are constrained in their 
value, either by collision geometry, by final state particle 
yields observed at RHIC, or/and by particle correlations. 
For example, the final yield of strangeness ds/dy and of 
entropy dS/dy, and thus s/S are known from an anal- 
ysis of particle production: as function of centrality at 
RHIC Q, and as function of reaction energy from top 
AGS, SPS to top RHIC energy J^. 

We will use these results in two ways. We compute, 
following the temporal evolution, the final state s/S ratio 
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which we expect to converge at RHIC to s/S ~ 0.033. 
We need to specify the initial value at time To for this 
variable, and this value is chosen to be compatible with 
the peripheral reactions. The entropy yield dS/dy, which 
we assume is conserved during the evolution of QGP, 
determines, for a known initial volume dV(r — To)/dy, 
the entropy density a = (dS/dy)/(dV/dy). We then can 
obtain, from standard properties of QGP fitted to the 
lattice results, the initial temperature Tq/T c ~ 3-4. This 
temperature decreases as volume expands with r given 
that the entropy is preserved. 

We will show that the two physical observables, s/S 
and 7s, we address are largely independent of the model 
dependent details of the initial conditions. Said dif- 
ferently, our important finding is that the two global 
strangeness observables s/S and "f s appear to penetrate 
back only to about 2 fm/c after the reaction has begun, 
and do not probe earlier conditions in the QGP phase. 
The physical reason for this is, of course, that once chem- 
ical equilibrium is approached, one looses much of the 
event memory with regard to intensive physical observ- 
ables. We will further see that, in cases we studied, that 
did not quite reach chemical equilibrium in the QGP 
phase, this is also true, i.e., there is little sensitivity to 
what exactly happened to light quarks and gluons in the 
first 2 fm/s. 

The reason for this is that there is a strong correlation 
between volume, temperature and degree of QGP (gluon 
and light quark) chemical equilibration as we already dis- 
cussed above. Repeating the argument differently, we can 
say that when fewer gluons at fixed entropy are in given 
volume, temperature has to be larger. Thus, any decrease 
in the production of strangeness in gluon fusion due to 
absence of gluons is compensated by the greater specific 
rate per colliding pair due to greater ambient tempera- 
ture. Hence, also when we do not quite reach chemical 
equilibrium in QGP, be it due to large impact parameter 
or low reaction energy (chemical noncquilibrium at lower 
energies is not explored in this paper), there is little if any 
dependence of thermal yield on initial conditions, and the 
results we arrive at regarding near chemical equilibration 
are extraordinarily robust. 

However, the thermal strangeness production does de- 
pend on the degree of initial state strangeness equilibra- 
tion, simply because, if the initial yields were chemically 
equilibrated to start with, there would be as much pro- 
duction as annihilation of strangeness and any temporal 
evolution is driven by the time dependence of the evolu- 
tion dynamics. We will take as a measure of the pre-QGP 
thermal phase strangeness production the specific per 
hadron multiplicity yield of strangeness observed in most 
peripheral RHIC reactions. This is typically s/S ~ 0.016 
at t = To . This choice allows to reproduce the observed 
value s/S = 0.019 attained in most peripheral nuclear 
reactions at RHIC [7(, with participant number about 
(A) = 6.3. Clearly, with s/S — > 0.033 in most central col- 
lisions, the implication of this choice is that the thermal 
process enhances total specific yield by factor 1.9 ± 0.3. 



As centrality and/or reaction energy decreases, there is a 
gradual decrease of this enhancement, and as the energy 
is increased (LHC) this enhancement rises somewhat. 



D. Strangeness production in thermal collisions 

We follow the established methods of evaluating ther- 
mal strangeness production (38| . expanding our earlier 
more schematic model (24| . However, considerable sim- 
plification arises since we focus attention on the specific 
yield of strangeness per entropy. 

In the local (comoving) frame of reference, the rate of 
change of strangeness is due to production and annihila- 
tion reactions only: 



1 ds 
Vdr 



1 ds 



88 



\&t) Mf?^ + Pq {t)pS){°)f 

-p s {t)p- s {t)(av) s r^. (3) 

The thermally average cross sections are: 

, v _ [ d 3 Pl J d 3 p 2 a 12 v 12 f(p 1 ,T)f(p 2 , T) 
( ™ rel)T " JcPp 1 Jd^ 2 f{p 1 ,T)f{f 2 ,T) ■ W 

f(pi,T) are the relativistic Boltzmarm/Juttner distribu- 
tions of two colliding particles i = 1, 2 of momentum p,-, 
characterized by local statistical parameters. 

A convenient way to address the dilution phenomena 
acting on the density of strangeness p s = s/V due to 
rapid expansion of the QGP phase, Eq. (J3j) , is to consider 
the proper time evolution of the specific strangeness per 
entropy yield: 



d s _V 1 ds 
d^S ~ ~S V dr" 



(5) 



The entropy S in a volume element is unchanged, as vol- 
ume grows and temperature drops: 



4tt 2 

S = V—g(T)T 3 = Const., 



(6) 



where we consider the quark and gluon degrees of free- 
dom along with their QCD corrections: 



2 S 8 C I 1 



15a. (T) 

47T 



--2 s 3 c n f 1 — 

4 V 2l7r 



(7) 



We use as the number of quark flavors rif ~ 2 + 
7 s 0.5z 2 K 2 {z), where z = m s /T. The terms proportional 
to chemical potentials are not shown in the expression 
for entropy, since p/ttT <C 1 at RHIC and LHC. 

We have used, here, the lowest order QCD correc- 
tions to the effective degeneracies, since these describe 
well the properties of QGP phase obtained on the lat- 
tice when the value of a s (T) used is as described 
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below, see Eq. (TT2")) . The agreement one sees for thermo- 
dynamic variables, such as E, P, S, with the lattice re- 
sults is very remarkable, including the temperature range 
near to the phase boundary. Thus, the use of constraint 
Eq. ^ to evaluate the time dependence of temperature, 
considering also that the third root of entropy is consid- 
ered, should yield precise enough results. We note, in 
passing, that we used as specified in Ref. [34[ the addi- 
tional terms A in entropy, arising from differentiation of 
the implicit temperature dependence of g, Eq. ((7|) enter- 
ing the partition function. 

In order to use the detailed balance which relates pro- 
duction and annihilation reactions, it is convenient to 
introduce the invariant rate per unit time and volume, 
A 12 ^ 34 , by incorporating the equilibrium densities into 
the thermally averaged cross sections: 



.4 



12^34 



1 



1 + A 



(8) 



$1,2 = 1 f° r the reacting particles being identical bosons, 
and otherwise, 6x2 = 0- Note also that the evolution for 
s and s in proper time of the comoving volume element 
is identical as both change in pairs. 

We find that the temporal evolution of s/S, in an ex- 
panding plasma, is governed by: 



d s 
~cfrS 



A99- 



(S/V) 
A q ^ 



+■ 



(S/V) 



bl(r)-l 2 s (r)] 
' K 2 (r)-7, 2 (r)] 



(9) 



When all 7^ — > 1, the Boltzmann collision term vanishes, 
and equilibrium has been reached. The value arrived at 
for the observable s/S depends on the history of how the 
system evolves and, eventually, reaches equilibrium. 

In order to be able to solve Eq. (9]), we need a rela- 
tion between s/S and 7 S . This is obtained combining 
strangeness density Eq. ([2]) and entropy Eq. ([6]) : 
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m a /T. 



(10) 



In the initial period, gluons and quarks have not reached 
chemical equilibrium, thus the actual numerical integrals 
of Bose and Fermi distributions of the type Eq. ([I]) , de- 
pendent on the values 7 g>g are employed instead, which 
modifies the result seen in Eq. (fl"0|) . 

The degeneracies we have considered in Eq. ([5]) for the 
entropy did include the effect of interactions, and thus, 
we have to allow for the interaction effect in the strange 
quark degeneracy as well: 



2 S 3 C 1 



ka s {T) 



(11) 



The value of k — 2 applies to massless strange quarks. 
At T = (or said differently, for m ^S> T) the early study 
of quark matter self-energy suggests that k — > [39| ■ We 



will present, in figure below, results for varying the 
value of k, and the reference value we use in our other 
studies is k = 1. We believe that this approach allows us 
to explore the general behavior of the interactions effect 
on the strangeness density, a more detailed study is not 
possible today. 



E. QCD parameters 

We evaluate A"^ ss and A q ^ ss employing the avail- 
able strength of the QCD coupling, and range of accepted 
strange quark masses. In our evaluation of strangeness 
production, in order to account for higher order effects in 
quark and gluon fusion reactions, we introduce a multi- 
plicative K — 1.7-factor. The known properties of QCD 
strongly constrain our results However, it turns out 
that the range of strange quark masses remains suffi- 
ciently wide to impact the results and we discuss this 
further in subsection IIIIDi see figure [7] below. We em- 
ploy when not otherwise stated the central value from 
a recent PDG evaluation, m s (fi — 2GeV) = 0.10 GeV 
which remains uncertain at the level of 25% at least. In 
fact, since our results were obtained, a more recent Parti- 
cle Data Group study recommends a 10% smaller central 
value of m s [40j |. 

We compute rate of reactions employing a running 
strange quark mass working in two loops, and using 
as the energy scale the CM-reaction energy fi ~ *Js. 
Since the running of mass involves a multiplicative fac- 
tor, the uncertainty in the mass value discussed above 
is the same for all values of [i. Some simplification is 
further achieved by taking at temperature T the value 
/i ~ 2ttT which is the preferred value of the thermal field 
theory, and agrees with the value of the reaction energy. 
This means that we use m s (T) = m s (fi = 2nT) with 
m s (T = 318 MeV) = 0.1 GeV. The actual temperature 
T(t) and thus time dependent values of the strange (and 
charm) quark mas will be always shown in the top panel 
of figures describing the evolution of the properties of the 
system considered. 

The strength of the QCD couping constant is today 
much better understood. We use as reference value 
a s (fi = rrizo) = 0.118, and evolve the value to applica- 
ble energy domain fi by using two loops. Our ability to 
use perturbative methods of QCD to describe strangeness 
production, a relatively soft process, derives from two cir- 
cumstances: 

a) the reaction processes which change yield of 
strangeness can compete with the fast v±_ > 0.5c ex- 
pansion of QGP only for T > 220 MeV, for lower tem- 
peratures the strange quark yields effectively do not 
change (strange quark chemical freeze-out temperature 
in QGP). Using the relation [i = 2irT, this implies that 
all strangeness yield evolution occurs for /1 > 1.4 GeV. 

b) Because of the magnitude a s (ji = m z o) = 0.118, one 
can quite well run a s to the scale of interest, /1 > 1.2GeV. 
As we see in figure [TJ this means that the strength of 
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fi [GeV] 

FIG. 1: The running QCD couping constant a s (fi) fixed to 
Q s (/i = m z o) = 0.118 (solid line) and several alternative 
strength scenarios excluded today by the experimental mea- 
surement (dashed lines). 



the interaction remains a s < 0.5. We also note that 
had the strength of a s (/j, = rn z o) been 15% greater, 
strangeness production could not be studied in pertur- 
bative approach. 

We next express a s (ji) (solid line in figure [T|) as func- 
tion of temperature by the conditions a s (T) = ct s (ji = 
2itT) This leads to the expression (see also section 14 
in Hi): 



u s (T c ) 



l + C\n{T/T c ) 



T < 6T C , 



(12) 



we consider, in the central rapidity domain: 



dV , , . dz 



(13) 



r— Const. 



dV/dy is the normalization factor for the particle yields 
we are measuring in an interval around central rapidity. 
The transverse expansion is described by the transverse 
size A± (r) . We further need to associate with the domain 
of rapidity dy a geometric region at the source dz, from 
which particles emerge. 

To accomplish this, we recall the space-time rapidity 
of the scaling Bj0rken hydrodynamical solution: 



y 



1. t + z 
— In . 

2 t-z 



(14) 



We see that y=0 corresponds to z = 0. In particular, if 
the transverse extend of the fireball is large, the Bj0rken 
space-time rapidity relation prevails. 

We need this relation not at a fixed laboratory time t 
but at some fixed proper time t: 



T = V< 2 - Z 2 . 

We eliminate in Eq. (Tl4|l t using Eq. (Tl5|) : 

z = rsmhy, — — rcoshy. 
dy 



(15) 



(16) 



A± is the transverse to scattering axes size of the evolving 
QGP. For nearly homogeneous expanding bulk matter 
one can assume: 



with C = 0.760±0.002, a s (T c ) = 0.50±0.04 at T c = 0.16 
GeV. We stress that Eq. (fT2"|) is a parametrization corre- 
sponding to the result shown in figure [TJ Only one loga- 
rithm needs to be used to describe the two loop running 
with sufficient precision, since the range we consider is 
rather limited, 0.9T C < T < 6T C . 



F. Expansion and cooling of QGP 

We separate, in our work, the issue of strangeness pro- 
duction from the even more complex and less understood 
questions about the time evolution of the QGP. We as- 
sume that there is some active volume at average temper- 
ature T, in which the strangeness is 'cooked'. We derive 
the time dependence of local temperature from the hy- 
pothesis of a conserved entropy content and a reasonable 
model describing the volume evolution in time. This is 
arrived at using a hydrodynamically inspired model. 

The volume at hadronization is an implicit observable. 
All particle yields at hadronization are normalized with 
a volume factor. Thus, our expansion model must be re- 
alistic enough so that the hadronization conditions are in 
agreement with data, and that the impact parameter de- 
pendence is reproduced. In the geometry inspired model 



A ± = ttRKt). 



(17) 



However, if the matter is predominantly concentrated 
near a narrow domain of width <i, we consider: 



A ± = n[Rl(r)-(Rl(r)-df] 



= 2nd 



(18) 



At central rapidity, we consider quantitatively the two 
evolution scenarios, denoted here-forth as models VI and 
V2. VI will be the most simple bulk homogeneous ex- 
pansion while V2 simulates a transverse donut, it corre- 
sponds to expansion with a cold hole of matter in fireball 
(axial) center: 



VI: ^f=-KRl{r)r, 
dy 



V2: ^ = 2nd 
dy 



with 



R ± (t) = R + v{r)dT 



(19) 
(20) 

(21) 
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Any model of transverse matter expansion dynamics 
v(t) is constrained by the transverse mass shape of pro- 
duced particle spectra, too large transverse velocities 
would produce too hard spectra. Accordingly, a hydro- 
inspired shape is assumed: 

2 

v(t) = V max - arctan[4(r - t )/t v ]. (22) 

Values of v max we consider are in the range of 0.5-0. 8c, 
the relaxation time t c ~ 0.5 fm, and the onset of trans- 
verse expansion tq was tried in range 0.1-1 fm. None of 
these parameters matters for what follows as long as one 
does not employ aberrant values. 

The initial size Rj_ is assumed, in what follows, to be 
R± = 5 fm for 5% most central collisions. When we study 
centrality dependence, we will show results for a series of 
centralities decreasing the transverse dimensions R± and 
d by factor Jr = 1.5 in each step. We further scale en- 
tropy value with fs = f^ 2 - This assures that the depen- 
dence of entropy on the participant number dS(A)/dy in 
the final state follows the relationship, 

^SL^-l), (23) 
dy 

obtained from the impact parameter dependent fit to the 
RHIC impact parameter results Q. 

An overview of the resulting volume dynamic behavior 
is given in figure O the top panel applies to RHIC with 
y / SNN = 200 GeV, the bottom panel presents a paral- 
lel study for LHC with the assumed four times greater 
entropy content. The solid lines are for transverse homo- 
geneous volume (VI model) expansion, and dashed lines 
correspond to a transverse region of thickness d = 3.5 fm 
(V2 model). 

Three different centralities were considered with R± = 
3, 5 and 7 fm. For the second model of transverse ex- 
pansion, the transverse size d is scaled with Rj_/5 fm. 
Thus, d = 2.1 fm for R ± = 3fm and d = 4.9 fm for 
R± = 7fm. Similarly, entropy content, assumed to be 
dS/dy = 5000 at RHIC and dS/dy = 20, 000 at LHC for 
R± = 5 fm, is scaled to values dS(R± — 3im)/dy = 1300 
and dS(R± = 7fm)/dy = 10,500, and correspondingly, 
4 times greater values for LHC. 

The temporal expansion of the volume is followed till 
T = 140 MeV is reached. In general, the maximum vol- 
ume at LHC is thus 4 times greater compared to RHIC. 
The expansion time is correspondingly longer, with RHIC 
taking 6.5fm to freeze-out for R± = 5fm, the LHC lifes- 
pan is lOfm. The QGP lifespan increase by as much as 
60% at LHC, when comparing to RHIC, if the assumed 
initial entropy production is indeed increased by factor 4. 
From the perspective of strangeness production, this is 
one of the more interesting changes comparing RHIC to 
LHC. 

Given the volume as function of r and the associated 
conserved entropy content, we can evaluate the prevailing 
temperature T for any given quark and gluon chemical 
yield condition 7g, s ,g- The solid lines in top panels of the 




0.1 1 10 

r [fm/c] 



FIG. 2: (Color online) QGP Volume related to central rapid- 
ity, dV/dy as function of proper time r. Top panel is for RHIC 
with reference entropy content dS/dy — 5, 000 (central lines), 
while bottom panel is for LHC with 4 times greater entropy 
content dS/dy = 20,000 (central lines). Three centralities 
are considered, with the middle thicker lines corresponding 
to R± — 5 fm and the upper/lower lines corresponding to 
R± = 7, and, respectively, R± = 3 fm/c. Solid lines are 
for VI model with transverse homogeneity, dashed lines for 
V2 model of a transverse shell with widths (top to bottom) 
d — 2.1, 3.5 and 4.9 fm. The volume expansion is shown in 
the figure up to T = 140 MeV. See text for more details. 



following figures [3] to [7] show this result, in the figures [3] 
and 21 on left for the VI model and on right for the V2 
(donut) model. The assumed 7 g is presented as dashed 
line. In some of the top panels, we also show by the 
dotted line the time dependence of the applied transverse 
velocity, v± , see Eq. (|2"21 . 



III. RESULTS ON STRANGENESS 
PRODUCTION 

A. The benchmark results for RHIC 

We present our results for strangeness production in 
QGP in Figs. [3H3 In Figs. [3J (RHIC) and H (LHC) 
we show the centrality dependence, with the two volume 
models, see Eqs. (fl"9l and l20l) . corresponding to columns, 
with VI on left and V2 (donut model) on right. In the 
following Figs. [SJ-Owe explore the dependence on the as- 
sumptions made. Here we show RHIC results on left and 
LHC results on right. In all Figs. [3]-[7]we show three pan- 
els above each other. As noted already, we show in the 
top panel, by solid linc(s), the model time-temperature 
profiles. The experimental observables are shown as solid 
line(s) in the middle panel (7^) and in the bottom panel 
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(s/S). The other lines illustrate as appropriate the key 
inputs used to obtain these results. When several lines 
of the same type are present, we are presenting the im- 
pact parameter dependence, scaling the size and entropy 
content as discussed above. In general, the temperature 
is followed down to a freeze-out at Tj = 0.14 GeV. 

In figure [31 we show what we believe is the best j s (r) 
(solid lines, middle panel) and s/S(t) (solid lines, bot- 
tom panel) for RHIC at 100+100 GeV at varying reaction 
centrality. The solid lines in the top panels show T(r) 
for these 6 different centralities, with the lowest tempera- 
tures seen for the least central collisions, all temperatures 
continue to T = 0.14 GeV. The slight increase in the ini- 
tial temperature with increasing centrality is result of the 
scaling of initial entropy, which accommodates the ob- 
served change in dS/dy\{ beyond participant scaling, see 
Eq. ([23]). For all centralities (and below also for RHIC) 
we assume the same initial s/S(tq) — 0.016. All lines 
shown begin at To = 1/4 fm, where the initial temper- 
atures range T G (0.55,0.6) GeV. For the VI model, 
the range of r spans the interval Tf = 2.2 fm (most pe- 
ripheral) to Tf — 6.5 fm (most central). In the donut 
expansion model V2, this range is from r/ = 3 to 8 fm. 

In the top panel, we also show the growth with r 
of the transverse expansion velocity (dotted line), and 
the strangeness pair energy threshold 2m r s using running 
strange quark mass (dashed lines). We note that the 
temperature drops below this threshold for the most pe- 
ripheral reactions considered already at r = 1 fm/c, and 
this occurs for the most central reactions at r = 2fm/c, 
for model VI, and respectively, 1.5 and 2.8 fm/c, for 
model V2. Thus, high strength thermal strangeness pro- 
duction life span varies by as much as factor 3, depending 
on centrality, and the expansion model. 

In the middle panel, we show dashed the rise of the 
gluon occupancy 7 g which we employed. The quark oc- 
cupancy 7 g is following the same functional temporal evo- 
lution starting with 2 /3 smaller initial value and evolving 
1.5 times slower. Because gluons dominate strangeness 
production in QGP, we do not show 7 g explicitly. The 
dotted lines show how the volume evolves toward its max- 
imum value at freeze-out. Each line is normalized to 
unity at freeze-out. The actual value of the volume can 
be read of the figure [21 given the value of r. 

To summarize the key results: we see a gradual in- 
crease of strangeness yield with centrality, reaching near 
strangeness QGP equilibrium for the most central colli- 
sions at RHIC. We have checked stability of this result 
against variation of model assumptions. More detailed 
discussion will be presented further below, see e.g. Fig. [6) 



B. Strangeness production predictions for LHC 

We performed a similar evaluation of strangeness pro- 
duction at LHC, see figure [|] which follows the same pat- 
tern as figure [3J There are three modifications which 
were introduced when we consider LHC: 



(i) To account for the greater reaction energy, as al- 
ready discussed, we increase the entropy dS/dy by factor 
4, which implies an assumed increase in rapidity den- 
sity of hadrons by a similar factor. We assume that 
in elementary parton interactions the relative strength 
of strangeness and non-strange hadron production is 
unchanged and thus, we keep the initial relative yield 
s/S = 0.016 constant. Given the entropy yield increase, 
we implicitly assumed an increase in initial strangeness 
yield by a factor 4 at LHC compared to RHIC. 

(ii) In order to accommodate the greater transverse ex- 
pansion pressure, we increased the maximum transverse 
flow velocity which can now attain v± = 0.80c (dotted 
line, top panel figure 3]). 

(iii) We further assume that thermalization time has 
dropped from t = 1/4 fm at RHIC to r = 1/10 fm 
at LHC. However, inspecting the slowly changing initial 
state evolution, in figure^ there would be little change in 
our results, were tq to remain unchanged between RHIC 
and LHC. At this early time, To = 1/10, the value of 
7s (to) at LHC is similar to the situation at RHIC, com- 
pare the beginning of the solid line in middle panel of 
figures [3] and [4j 

This is so, since the magnitude of the phase space scales 
with T 3 and the initial temperature T(tq) is considerably 
greater at LHC: in the top panel, in figure [H we see 
that it reaches up to T = 1.25 GeV. For this reason, we 
have to show in the top panel (dashed lines) bm r s rather 
than 2m r s , in order to fit it visibly into the top panel of 
the figure, and this is the only difference in the display 
of LHC results, in figure [U as compared to the RHIC 
results, figure O 

We note that despite a much greater expansion ve- 
locity, the evolution time at LHC is significantly longer, 
with the most central collisions taking up to 30% longer 
to reach the freeze-out temperature, Tf = 0.14 GeV. The 
reader who prefers earlier freeze-out, at, e.g., Tf = 0.17 
GeV can evaluate the changes required by consulting the 
temperature profiles shown in the top panel. 

The different centralities, at LHC, are considered with 
the same scaling of the transverse size and entropy con- 
tent as we did for the case of RHIC. Comparing the 
left to right set of results (i.e., volume VI to donut V2 
expansion models) for LHC, we see a more developed 
chemical equilibrium of strangeness with clear evidence 
of (over)saturation of yields for a few centralities, see the 
bottom panel on right in figure [U There is greater final 
specific strangeness content at LHC than at RHIC, with 
visibly greater thermal production leading to strangeness 
(over)saturation. 

At RHIC, the thermal production raises the value of 
s/S from 0.016 to 0.028 for most central collisions (VI 
model of bulk volume expansion), while at the LHC the 
thermal production raises s/S from 0.016 to 0.032. We 
will discuss below what this increase means for the K/tt 
and other particle ratios. The same relative increase in 
s/S is seen in the model V2 of the expansion compar- 
ing RHIC and LHC. However, if the homogeneous bulk 



10 



0.03 - 



0.02 



i — mi 1 1 — i — i..t i i 




"1 u fTlf \ ' \ \ — I — I I I I 




10 



[fm/c] 



i — in 1 — i — i — i — i 




[fm/c] 



FIG. 3: (Color online) RHIC results. Top panel: solid lines: temperature T; dashed lines: running mass m r a {T)\ dotted line: 
the assumed profile of transverse expansion velocity v±(t). Different lines correspond to different centralities. Middle panel: 
Solid line(s) 7 S , which nearly coincide for different centralities; dashed, the assumed Js(t), dotted the profile of the assumed 
volume, \dV{r) / dy]/[dV (jf) / dy\ normalized by the freeze-out value. R±(tq) stepped down for each line by factor 1.5. The 
end points, at maximum r, allow to identify corresponding s/S for different centrality in the bottom panel. Right and left: 
Comparison of the two transverse expansion models, see Eqs. (| 19I20[) : left bulk expansion (model VI), right donut expansion 
(model V2). 



expansion applies at RHIC, but a donut type expansion 
arises at LHC, the increases in strangeness yield and lifes- 
pan would be more spectacular. Depending on its ex- 
pansion dynamics, LHC clearly harbors the potential to 
surprise us. 



Study of the dependence on initial 
thermalization condition 



An important question is how the value of the unknown 
initial conditions impacts the results we presented above. 
We have studied this question in depth in many different 
model approaches. The answer 'practically no depen- 
dence' is best illustrated in the figure where we show 
the more conservative volume expansion model VI results 
on left for RHIC, and and on right for LHC. We explore 
a wide range of initial gluon (and quark) occupancy 7 g , 
which for consistency with other figures is shown in the 
middle panel by dashed lines, the initial values we con- 
sider for glue occupancy vary as 0.1 < 7 g (To) < 2.1 in 



step of 0.5. The second of these lines, from the bottom, 
is the reference behavior we use in figures [3] and 2] 

We recall that, with 7 g (ro), we also vary 7 9 (to), which 
following the same functional temporal evolution start- 
ing with 2/3 smaller initial value and evolving 1.5 times 
slower. Note also that the scale, in top panel, varies be- 
tween RHIC (left) and LHC (right) cases the dashed lines 
denote 2m s on left, and 5m s on right. Since the initial 
value of s/S = 0.016 and dS/dy = 5, 000 on left for RHIC 
and, respectively s/S = 0.016 and dS/dy = 20,000 on 
right for LHC is set, there is a corresponding variation in 
T (top panel, left end of solid lines) and 7 S (left end of 
solid lines in middle panel). The final results for 7 s (t/) 
(right end of solid lines in middle panel) and s(rf)/S 
(bottom panel) are impressively insensitive to this rather 
exorbitant diversity of initial conditions at fixed entropy 
content. The spread in s/S(t) we see in the bottom panel 
could be seen as a wide line width. 

We conclude that strangeness cannot probe the very 
initial QGP conditions near tq — the memory of the ini- 
tial history of the reaction is lost, the system is opaque 
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FIG. 4: (color online) Case of LHC, see legend in figure [3] See text for discussion of differences with RHIC. 



for t < 2-3fm/c to the strangeness signature. On the 
other hand, and most importantly, for the study here un- 
dertaken, this also means that experimental observables 
are characteristic of the properties of nearly chemically 
equilibrated QGP. 

However, there remains dependence on the history of 
the QGP fireball in models in which gluon (and light 
quark) chemical equilibrium is not attained even in cen- 
tral reactions at RHIC [25l . [26| . Thus only if gluons 
in QGP did not approach the chemical equilibrium at 
Tf, a signature of this condition would be seen in the 
strangeness yield as is seen considering the right hand 
side of Eq. @. Correspondingly smaller value of s/S are 
then expected. 

Yet, our here presented results agree well with cur- 
rent RHIC strange hadron production results regarding 
the total strangeness yield. This constitutes indirect 
evidence for the achievement of light quark and gluon 
chemical equilibrium in QGP formed in the most cen- 
tral, highest A and highest energy RHIC reactions. On 
the other hand, the relatively small size and small lifes- 
pan of a QGP potentially also formed in peripheral col- 
lisions could hinder the achievement of QGP chemical 
equilibrium at hadronization. 



D. Fundamental uncertainties in strangeness 
production 

There are two QCD related uncertainties in our 
strangeness production study which we explore in turn 
in figures [5] and [7] 

1. the effect of interaction on the number of strange 
quark degrees of freedom, Eq. (fTTjl , 

2. the value of strange quark mass. 

We will now show that these uncertainties lead to ob- 
servable effects, in particular regarding the final value of 
s/S, and to a lesser extend also regarding the final value 
of 7 S (chemical yield equilibration). 

We first recall that 7 S is introduced in Eq. @ in order 
to relate the prevailing strangeness density to the chem- 
ically equilibrated density at temperature T. Given a 
strangeness yield, the value of 7 S depends on the prop- 
erties of QCD by the way of what the chemically equi- 
librated density is. The actual value of ~f s enters deci- 
sively into the kinetic equation of strangeness produc- 
tion. Eq. © shows that the smaller the value of 7 S , the 
bigger is the obtained change in the value of s/S. 

Since, in our approach, we used Boltzmann statis- 
tics for the strangeness degree of freedom, the effect of 
Pauli blocking on strangeness production is not consid- 
ered and in fact is a minor effect. The direct depen- 
dence of physical observables on j s arises from the pro- 
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FIG. 5: (color online) Model VI (volume expansion) at RHIC (left) and LHC (right) for 5% most central collisions. s/S 
(bottom panel) and 7 S (solid lines middle panel) as function of r, for widely varying initial gluon conditions (T in top panel, 
7g dashed middle panel), constrained to same entropy content. 



cess of strangeness reannihilation into gluons. We note 
that Eq. © , describing the change in strange quark yield, 
can also be written in the form: 



drS 



a s (pipT 2 -pip^ A )+a q (py s 



(24) 



where instead of 7i the actual densities of particles ap- 
pear. We see that for each particle both the equilibrium, 
and transient, density must enter in order for the system 
to attain chemical equilibrium as function of time. Con- 
sequently, the value of 7 S and hence the QCD correction 
in Eq. ((2j) matters. 

The figure O which follows the pattern of figure 
El illustrates this effect when the effective degeneracy, 
see Eq. ([IT]) , varies. We consider for values of k = 
2, 1.5, 1., 0.5 and 0. k = 2 is the perturbative effect 
seen for massless quarks, k — corresponds to no effect 
of interaction, when m s ^S> T. When k = 0, the strange 
quark degeneracy is largest, thus for a given strangeness 
yield s the value of "f s is smallest, and the production 
of strangeness biggest. Consequently, this value corre- 
sponds to the smallest 7 S in the middle panels of figure 
[5] and greatest value of s/S in the bottom panels. Other 
lines follow, and the middle solid lines (red on line ) cor- 
respond to k = 1, which we used to obtain the reference 
results presented earlier. 



A similar effect arises considering variation of strange 
quark mass as shown by dashed lines in top panels of 
figure [71 We vary by factor 2 the strange quark mass as 
is seen in the top panel of the figure \7\ on left for RHIC 
and on right for LHC. When the strange quark mass is in- 
creased, the equilibrium strangeness density is decreased, 
and thus, for a given strangeness yield s/S, the value of 
7 S is increased, which in turn reduces strangeness pro- 
duction strength. Smallest mass considered (dotted line 
in top panels of figure [7]) corresponds to the largest final 
value of s/S shown by dotted line in the bottom pan- 
els. What is more surprising is that the effect of mass 
variation cancels in the actual computed 7 S shown in the 
middle panels. The reason for this accidental cancella- 
tion is that for a larger mass the smaller value of final 
s/S solves for the same value of %, see Eq. (fTU|) . 



Since both m s , and the interaction effect on the strange 
quark effective degeneracy, see Eq. (flTf , are today not 
understood at sufficiently precise level, the appearance of 
a possible range of values at freeze-out, for both 7 S and 
s/S, in figures [5]and [3 signals (correlated) uncertainty in 
the understanding of the results at RHIC and predictive 
power for LHC. 



13 




FIG. 6: (color online) Study of of s/S evolution with r, on left for RHIC, on right for LHC, VI volume expansion model. 
Figure structure same as figure [5] Middle panel, solid lines: computed evolution of 7 S in the deconfined phase for values of 
k — 2, 1.5, 1, 0.5 and 0, see Eq. (|11[) ; Bottom panel: corresponding evolution of s/S. The lowest 7 S line corresponds to k — 
and the largest to k = 2, the opposite applies for s/S lines. 



IV. CONSEQUENCES FOR HADRON YIELDS 
AND THEIR EVOLUTION FROM RHIC TO LHC 

A. Strangeness and entropy 

The final value of s/S is the key result which prac- 
tically alone determines relative strange particle yields. 
It depends on the value of dS/dy which we start with. 
The value we take, dS/dy = 20,000 at LHC (4 times 
RHIC), is a guess arising from extrapolation of the en- 
ergy dependence of particle production. However, it fixes 
in effect the initial conditions, and leads to the range of 
values for s/S we obtained. We had found that there is 
continued growth of strangeness yield beyond the RHIC 
energy range where dS/dy = 5, 000, which suggests a fur- 
ther strong increase in the strange hadron yield we will 
address below. 

Variation of dS/dy amounts to variation of the pro- 
duced hadron yield. We expect that dS/dy cx dh/dy. To 
quantify this, we present, in figure [8l how the yield of 
charged hadrons relates to the entropy yield. These re- 
sults were obtained both for RHIC (left region, red on 
line) and LHC (right region) by the methods described 
in Ref . [§| , employing SHARE suite of programs [4l[ . In 



this evaluation we have available for each s/S (at RHIC 
and LHC) the corresponding strangeness occupancy j s . 
We further can fix the ratio of charge to baryon number 
as in the incoming nuclei, q/b = 0.39. 

A further difference between RHIC and LHC is that 
we take the thermal energy per baryon at RHIC to be 
E/b = 39.3 GeV and following Ref. [g E/b = 412 GeV at 
LHC. The consistent ranges for s/S at RHIC are 0.018 < 
s/S < 0.03 and, at LHC, we can address 0.018 < s/S < 
0.037, we cannot otherwise find a smooth match of QGP 
to HG-phase space within the physical range of phase 
space occupancies 7 S , jq. We note that this LHC choice, 
E/b = 412 GeV, limits the range of possible entropy yield 
to just below the range we explore in our present work, 
dS/dy = 20, 000, four times the RHIC range. 

We next explore the resulting final s/S and j s . We 
show these quantities, on left in figure [9l as function of 
the dS/dy input, and as function of the charged hadron 
multiplicity d(h~ +h + )/dy, on right. The two expansion 
models we consider are as before VI homogeneous ex- 
pansion (thin solid, lower line) and V2 donut expansion 
model (thick solid, upper line). We see a gradual rise of 
both s/S and 7 S as function of dS/dy which begins to 
saturate for dS/dy > 20, 000 , but at a rather high val- 
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FIG. 7: (color online) Study of s/S evolution with r for different m s , on left for RHIC, on right for LHC, VI volume expansion 
model. Figure structure same as figure [S] and [ll] Top panel: 5 (running) strange quark mass m r a . The middle of 5 lines being 
our standard reference value, see text for more detail. 



3000 




10- 4 dS/dy 

FIG. 8: (color online) The yield of charged hadrons d{h~ + 
h + )/dy for different values of dS/dy, left domain for RHIC 
and right upper domain for LHC. Solid lines: before weak 
decays, dashed lines: after all weak decays. 



ues beyond chemical equilibrium: the expansion is so fast 
that there is no time to reannihilate the very abundant 
strangeness before freeze-out. 

In Fig [9] on the right hand side, we see the RHIC 
(smaller dh/dy range) and LHC domains. Dashed lines 
are the charged hadron multiplicity after weak decays, 
in particular of neutral strange hadrons such as Kg. The 



negligible impact of charged Kaon weak decays to hadron 
multiplicity is not accounted for as this effect is small 
and experiment dependent considering the quasi-stability 
of K ± . In the RHIC domain, we have a nearly lin- 
ear rise of s/S (top panel) and 7 S (bottom panel) with 
d(h + + h~)/dy. In the LHC domain chemical (over) sat- 
uration of these quantities is clearly visible. The left, 
and right hand side of figure [H] have at large dS/dy and 
d{h + -\-h~) / dy, respectively, a similar appearance indicat- 
ing chemical (over)saturation with the increasing hadron 
yield. We recall that while the functional form of the fig- 
ure is correct, the normalization of the abscissa is some- 
what arbitrary for the case of LHC, which is based on an 
extrapolation of SPS and RHIC results. 

We next study the growth of strange hadron yield with 
s/S. In order to not confuse this with the issue of final 
hadronization volume of QGP, we will study particle ra- 
tios. The method of evaluating particle yields is, as just 
discussed in study of the charged hadron yield: SHARE 
program is asked to match micro-canonical conditions in 
the QGP phase to those of HG phase space. Given the 
statistical parameters, we can compute all particle yields 
within the statistical hadronization approach. 

We first consider, in the two top panels of figurefTUl the 
ratios A/pi + and K + /tt + as function of s/S, the attained 



15 



0.04 



0.5 



10 3 



n i i i i i 



n i i i i i i I 




\ 1 I I I I II 



_i I I 



_i I i 



10 4 

dS/dy 



10 5 



0.04 



m 
m 



0.03 



0.75 



i i i i i i i I i i i i i i i i i I i i i i i i i i i 




I I I I I I I I I I I I I I I I I I I I I I I I I I I - 








1000 2000 
d(h + + h-)/dy 



3000 



FIG. 9: (color online) s/S (top panel) and 7 S (bottom panel) for the most central 5% collisions as function of dS/dy (on left) 
and as function of d(h + + h~)/dy on right. The bulk transverse expansion model VI is the thin solid line (blue), the thick line 
(red) is the donut expansion model V2 with d = 3.5 fm. On right, solid lines: before weak decays, dashed lines: after weak 
decays (excluding K see text). We recognize the RHIC domain results by the smaller dh/dy range of results presented. 



specific strangeness per entropy at QGP freeze-out. Solid 
lines are before weak decays, and dashed lines present the 
corresponding results after weak decays. These enhance 
the pion yield, but as we see, even more the A yield. The 
K + yield is practically not changed: considering the long 
life span of the charged kaon it is common to present 
experimental results corrected for any decay. Thus, 
are considered, in our study, as if they were stable parti- 
cles. The RHIC results are the thicker lines and the LHC 
results are the thinner lines. 

The RHIC and LHC results for K+/7T+, 0/K+ and 
Sl~/K + are practically overlapping, since the influence 
of the difference in baryochemical potential is not mate- 
rial for this ratio. However, the greater final s/S yield 
expected in most central LHC collisions implies an in- 
creased yield at LHC compared to RHIC. The dotted 
lines in the second from top panel, in figure [TUl guide the 
eyes both in the RHIC domain, and at the higher s/S, 
LHC domain. 

We find, in this perhaps easiest to measure K + /7r + ra- 
tio, a noticeable and measurable relative yield increase. 
This prediction is important, since after dropping with 
energy at SPS, this ratio remained of the same constant 
magnitude (within error) at the much higher RHIC en- 
ergy domain. There is less change expected between 
RHIC and LHC for the A/pi + ratio, since the decreasing 
baryochemical potential is compensating, to a large ex- 
tend, the effect of the increase in strangeness yield. We 
note that there will be the opposite effect in the A/ir + 
yield. We see that the ratios 4>/K + and fi/K + are most 
sensitive to the change in s/S occurring between RHIC 
and LHC. The change in baryochemical potential dimin- 
ishes this sensitivity considerably in the S/K + ratio. Ex- 



cept for A/K+ the weak decays have a negligible impact 
on all ratios with K + considered here. 



B. Strange particle yields as function of centrality 

We have obtained the growth of strangeness yield in 
QGP with centrality and this allows us to explore partic- 
ipant dependence of strange hadron yields. Since we do 
not know well how E/b and E/TS changes as centrality 
changes, we choose to follow a different approach in order 
to present our results. We assume that the hadronization 
occurs at a fixed condition, and first we consider T = 140 
MeV and 7 9 = 1.6, corresponding to supercooled sudden 
hadronization. In a second step, we compare these re- 
sults to those obtained with T = 160 MeV and j q = 1, 
which is the equilibrium hadron phase space. The fixed 
value of statistical hadronization parameters as function 
of centrality follows the pattern seen in the analysis of 
RHIC results Q. 

To obtain the results, we follow a similar computa- 
tional scheme as developed in the study of s/S depen- 
dence, except that we now fix T MeV and 7 g , and not, 
e.g., E/TS. These conditions lead to similar hadroniza- 
tion conditions, but E/TS had some variability for small 
participant number, see figure 4 in Ref. [tJ- 

In the figure fTTT we show the results, on left for RHIC, 
and on right LHC. We follow a similar display scheme 
as in figure [TOl and show A/pi + , K + /tt + , followed by 
0/K+, A/K+, 5-/K+, fi-/K+, as function of partic- 
ipant number. The thinner (blue) lines are for the VI 
model of volume expansion, and the thicker (red) lines 
for the V2 model of donut expansion. Dashed lines in- 
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FIG. 10: (color online) Relative particle yields as func- 
tion of s/S: from top to bottom A/ir + , K + /tt + , and 
0/K+, A/K+, E~/K + , fT/K + . Solid lines (blue) are the 
primary relative yields, dashed lines (red) give the yields after 
weak decays (K + is not decayed). Thick line with s/S < 0.03 
are for RHIC, and thin lines with s/S < 0.037 are for LHC 
physics environment. Dotted lines guide the eye for the RHIC 
K+/7T+ ratio. 



elude the weak decays, which increase the A/pi + ratio 
and decrease the K + /ir + ratio. 

We see, in figure [Til a steady, but rather slow, in- 
crease with centrality of the relative strange hadron, to 
non strange pion yield. The decay pions tend to flatten 
the K + /ir + ratio at RHIC. We obtain predict more rise 
of (j>/K + than is observed at RHIC [42|, but our variation 



is within the error bar of the experimental, nearly con- 
stant result, <fi/K + = 0.15±0.03. As discussed earlier, the 
overall increase in s/S expected at LHC (right side) com- 
pared to RHIC (left side) explains the noticeably greater 
relative K + /tt + ratios at LHC. The most noticeable rise, 
with centrality, is expected when the ratio has the largest 
disparity in strangeness content. There is no central- 
ity dependence expected when there is no difference in 
strangeness content, such as in A/K + , where the weak 
decay produces very small centrality dependence. 

One of the interesting questions is how sensitive is the 
study of these particle ratios to the hadronization condi- 
tions. In figure we compare the sudden hadronization 
at T = 140 MeV and j q = 1.6 with the hadron phase 
space equilibrium model at T = 160 MeV and j q — 1 
forming the ratio of particle yield ratios, that is the re- 
sults obtained for sudden hadronization are divided by 
those obtained for the equilibrium HG phase space. The 
panels follow the same particle ratios as in figure QTJ on 
left we show RHIC and on right LHC results. 

The value of "f^(T = 160,7^ = 1) varies as func- 
tion of centrality 1.26 < 7^ < 0.88. In general, the 
higher hadronization temperature assumed in this chem- 
ical equilibrium case favors the yield of the more massive 
hadron. Thus, ratios of heavy to lighter particle evalu- 
ated at the same value of s/S as shown in figure [TT1 are 
bigger for the 7J 1 = 1 since we took a higher T = 160 
MeV chemical equilibrium freeze-out value. 

We note that the nonequilibrium model is much bet- 
ter explaining available multistrange hadron data. The 
equilibrium model needs to be amended to explain the 
enhanced fl yield. The way out from this dilemma if 
one insists on HG equilibrium could be a multi-freeze- 
out temperature interpretation. The freeze-out tempera- 
ture of the Q~ in the equilibrium freeze-out model would 
need to be noticeably higher than that of the bulk of 
strange particles. We note that such multi-freeze-out 
models could experience systematic difficulties as func- 
tion of centrality. Moreover, considering figure [T2l in the 
chemical equilibrium interpretation of hadron production 
a separate freeze-out for both £l~ and S~ will be required 
at LHC. Aside of being unpalatable, such a multi-freeze- 
out HG equilibrium model will not describe fluctuations 

well m. 



C. Thermal charm at RHIC and LHC 

The direct initial high energy parton collisions domi- 
nate production of the massive flavor, charm and bottom, 
and a thermal process seems to be of no interest. How- 
ever, there are two questions we can investigate: 

a) considering that 7 C in the deconfined phase can be as 
large as 7 C ~ 100 is there any significant thermal annihi- 
lation of charm in the QGP evolution? 

b) how large is the thermally produced charm yield and 
can it lead to chemical equilibrium of charmed quarks? 
The interest in question a) needs no further discussion. 
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FIG. 11: (color online) Relative strange particle yields as function of participant number A, left RHIC and right LHC. From 

top to bottom, A/tt + , K + /tt + , and <j>/K + , A/K+, H"/K + , JT/K+. Solid lines primary relative yields, dashed lines (relative) 

yields after all weak decays (not shown when difference is within line width). Thin lines (blue), model VI (volume expansion) 

and thick lines (red), model V2 (donut expansion). Results are for supercooled sudden hadronization. 



Question b) is of interest since the directly produced 
charmed quarks are in principle not easily thermalized, 
while the thermally produced charmed quarks emerge 
naturally in a momentum distribution, imaging their 
'parent' particles thermal momentum distribution. Con- 
sequently, the thermally produced charmed quarks pro- 
vide a solid thermal lower limit for the yield of charm, 
with the directly produced charm contributing to thermal 
distribution after charm has been subject to collisions re- 



quired for thermalization. 

We see the results of this study in figure [T31 on left 
for RHIC, on right for LHC. The top panels, as usual, 
presents the temperature and, charmed quark mass, 
scaled with factor 1/10 at RHIC and 1/2 at LHC (on 
right). The middle panel presents 7 g , the charm phase 
space occupancy, j c and j c obtained solely by thermal 
processes. Similarly, in the bottom panels, the dashed 
lines are the thermal yield at RHIC and LHC, while the 
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FIG. 12: (color online) Ratio of relative yields obtained for sudden and equilibrium hadronization (see text) as function of 
centrality. Left RHIC, and right LHC. From top to bottom, A/tt+, K+/tt + and 4>/K + , A/K+, H~/K+, fT /K+ . Solid lines 
primary relative yields, dashed lines (relative) yields after all weak decays (when absent, no difference within line width with 
solid lines). Thin lines (blue), model VI (volume expansion) and thick lines (red), model V2 (donut expansion). 



horizontal lines are the (little) evolving c/S yields in- 
cluding the directly produced charm. The direct charm 
production, at RHIC, is expected to be 600 times greater 
than the thermal process. At LHC, the higher initial 
temperature, but unchanged specific direct yield c/S, in 
parton collisions, suggests that thermal production is just 
factor 90 times smaller. However, this factor depends 
on good understanding of both processes and the initial 
conditions and surely cannot be fully trusted. Moreover, 
there is the possibility that the directly produced charm 
at LHC may not well thermalized, in which case an ap- 
preciable fraction of thermal charm yield could indeed 



originate in thermal reactions. 

The answers to the opening questions thus are 

i) There is no visible charm reannihilation, see the nearly 
horizontal lines at c/S = 3 • 10~ 3 (left, RHIC) and c/S = 
6 • 10~ 3 (right, LHC), in the bottom panel; note that the 
direct charm yield, we assumed implicitly, is obtained by 
multiplying the c/S yields with dS/dy = 5, 000 on left 
for RHIC, yielding dc/dy\xmic = 15, and with dS/dy = 
20, 000 on right for LHC, yielding dc/dy\ h HC = 120. 

ii) In the middle panel, we see that thermal production 
alone (dashed lines) oversaturates the charm phase space, 
the thermal produced charm phase space occupancy 7* 11 
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FIG. 13: (color online) Left RHIC and right LHC for charm production. Figure structure same as figures[3H7] Top panel: solid 
lines T, dashed lines, running m r c , scaled with 10 for RHIC on left, and with 2 on right for LHC. Middle panel: dotted line 7 g , 
solid lines the computed total charm y c , dashed lines 7 C corresponding to thermal charm production. Bottom panel: specific 
charm yield per entropy, solid lines for all charm, and dashed lines for thermally produced charm. 



(dashed lines middle panels) cross the gluon dotted j g 
line at around 4 fm/c for RHIC (corresponding to T ~ 
0.175 GeV) and at 5 fm/c (corresponding to T ~ 0.20 
GeV) at LHC. 



V. SUMMARY AND CONCLUSIONS 

We have studied the thermal QGP based strangeness 
production at RHIC and LHC, and have interpreted the 
observed final s-yield in terms of our theoretical knowl- 
edge about the properties of the QGP phase. Our aim 
has been to understand how the overall final state strange 
quark flavor has been produced, and to study in detail 
the mechanisms behind strangeness enhancement. As a 
further objective we have explored the impact of high 
strangeness on the strange hadron yields. 

Our results suggest that strangeness enhancement 
could be studied considering: 

ps _ S/S | ccn t r al _ s/S(r f ) 

n CP = ~Tc\ — To? — \' ( °) 

S/O | peripheral S/b(T ) 

The central strangeness yields is just the final value we 
find at freeze-out, combining the initial yield with the 



additional thermal production. The peripheral yield is 
the initial value before thermal strangeness production 
begins. Our study shows that R^ P € [1.6,2.2], with the 
precise result depending on details such as strange quark 
mass, see figure [7j reaction energy and dynamics of ex- 
pansion, see figure [2J 

More generally, to separately consider s and S we can 
use d(h + + h~)/dy as a measure of entropy dS/dy con- 
tent, see figure [8] Instead of the total strangeness, on 
may consider enhancement of individual (multi)strange 
particles, which we discussed in depth both as function 
of the achieved s/S, in subsection IIV Al and as function 
of ccntrality A at fixed hadronization condition for all 
centralities, in subsection IIV Bl 

We note that the overall growth of the enhancement 
of the strangeness yield with centrality, at the level a 
factor 1.6-2.2 is accompanied by a further enhancement 
of multistrange hadron yields, as is seen comparing the 
yields of multistrange hadrons with the yield of kaons. 
One can show that the </>/K + ratio is mainly dependent 
on value of s/S and not on hadronization temperature 
(see appendix B2i 43]) when the strangeness conserva- 
tion constraint is implemented. This effect is unique to 
<fi and arises since ~ 2mK- 

Aside of centrality dependence, we have explored, 
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within the framework of our model, the extrapolation 
from RHIC to LHC physics environments. More gener- 
ally, we have presented, in figure [H the reaction energy 
dependence by considering the central rapidity s/ S yields 
as function of dS/dy. 

One of the interesting results obtained is the approach 
to chemical strangeness equilibrium in the deconfincd 
QGP phase formed in most central and high energetic 
RHIC reactions. The evidence for this is implicit in the 
experimentally reported yields of strange hadrons, which 
lead to values of specific strangeness per entropy at the 
level of s/S ~ 0.028 % Our study of QGP based kinetic 
strangeness production provide an explanation of this re- 
sult, both, the value of s/S and as function of centrality. 

Our present study further shows that the proximity of 
chemical strangeness yield equilibration in QGP formed 
at RHIC and LHC, and the effective opacity of QGP to 
this signature, is the reason that considerably less so- 
phisticated models of QGP evolution which we, and oth- 
ers, have considered are equally successful in the study 
the strangeness production, as long as these models yield 
conditions near to chemically equilibrated QGP. 

Given the near chemical equilibration at RHIC, 
and within the models considered, we obtain some 
strangeness over saturation at LHC. Moreover, for the 
most central 5% reactions there is no relevant dependence 
of strangeness production on initial conditions prevailing 
in the reaction. We have demonstrated this in a picture- 
book fashion, see figure [5l where, for a wide range of 
initial conditions, the same final strangeness yield and 
equilibrium condition arises after At =2-3 fm/c. 

On the other hand, the more peripheral reactions do 
not saturate the phase space, in that both 7^ GP < 1, 
and s/S < 0.03. Thus, the peripheral yields, being sensi- 
tive to the initial conditions, allow exploration of physical 
conditions in the QGP prior to the onset of chemical re- 
actions. Therefore, our results for most peripheral reac- 
tions are also somewhat dependent on model assumptions 
about initial state and evolution dynamics. 

We have shows, in figures [3] and 2J the impact parame- 
ter dependence that arises in two volume expansion mod- 
els. Since the analysis results presented in Ref. 0] have 
been used to fine tune the dynamical evolution model at 
RHIC, there is good qualitative agreement with experi- 
mental results. Our objective has been here to learn how 
to extrapolate the dynamics of strangeness production to 
the LHC environment. The gradual rise of strangeness 
yield s/S with dS/dy, seen in figure O is reminiscent of 
the rise of s/S with reaction energy obtained in a analysis 
of particle yields obtained at SPS and RHIC at different 
reaction energies [37j • 

We have made a (conservative) prediction regarding 
the increase in the K + /n + ratio at LHC compared to 
RHIC, sec figure [10] Though an important result is that 
we expect an increase at LHC, we note that even a greater 
increase is possible, signaling even greater values of s/S, 
depending on both: 

1) the dynamics of the volume expansion — this can 



enhance the strangeness oversaturation of the final QGP 
state, see figures [3] and 0] 

2) QCD details, such as strange quark density including 
QCD interactions, and (still not well understood) strange 
quark mass, see figures [5] and [7] 

In any case, we believe that this simple observable will 
show again strangeness production growing faster than 
entropy production, its increase is directly coupled to an 
increase in s/S. We note again that in our study the 
increase of s/S with centrality implies that <^>/K + also 
increases with centrality. 

A natural result is the finding of the chemical yield 
equilibration of strangeness in the QGP formed in the 
5% most central top RHIC energy reactions. This leads 
to a better understanding of the resulting oversaturation 
of the hadron phase space by strangeness. The mag- 
nitude of this effect, dependent on the temperature of 
hadronization, can be considerable. This can be eas- 
ily seen considering the magnitude of s/S in both QGP 
and HG phases. The final state hadrons formed far-off 
chemical equilibrium cannot significantly adjust chem- 
ical composition, considering the rapid breakup of the 
fireball, during the period of about 1-2 fm/c prior to on- 
set of the free flow. Thus, our finding is that strangeness 
rich QGP enhances decisively the yields of multistrange 
hadrons. This phenomenon is more accentuated consid- 
ering charmed hadrons containing strangeness, a topic 
under current investigation. 

Furthermore, using the methods developed here, we 
have considered thermal charm production. At LHC, we 
find a nearly physically relevant thermal charm produc- 
tion, but not at RHIC. However, the thermal process we 
consider is able to produce enough charm to oversaturate 
the final state at both RHIC and LHC, see figure[T3l This 
also shows that the direct parton collision based produc- 
tion at RHIC leads to extraordinarily large values of j c . 
The chemical nonequilibrium of charm is thus more pro- 
nounced than that of strangeness. 

In conclusion: The totality of our results shows that, 
as function of entropy yield dS/ dy (equivalently, the re- 
action energy of A\-Ai collision) and geometric reaction 
size (impact parameter dependence, participant number 
A) , the phenomenon of strangeness enhancement is well 
described by the mechanism of QGP based thermal gluon 
fusion strangeness production. We find both, as function 
of centrality, and energy, a gradual increase in specific 
strangeness yield, which agrees with all available experi- 
mental results. We find that, as function of energy, this 
continues from RHIC to LHC increasing our hopes for a 
more clear strangeness signature of deconfinement. 
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